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The resonance model is used to analyze the uj- and 0-meson angular distributions in proton-proton 
collisions at y/s — 2.83 and 2.98 GeV. The assumption of dominant contributions from iV*(1720)§ + 

and iV*(1900)§ + resonances which both have, according to the ttN scattering multichannel partial- 
wave analysis and/or quark models predictions, dominant p 1( / 2 Nuj decay modes yields the right 
pattern of the ui angular distribution at y^s = 2.83 GeV. The angular distribution at y^s = 2.98 GeV 
can be reproduced assuming the dominance of iV* (2000) § and iV* (1900) § . The experimental 0- 
meson angular distributions do not shown any asymmetry which requires the existence of a massive 
' negative-parity spin-half resonance. This resonance could be identified with the iV* (2090) ^ ■ 

(N ■ 

■ PACS numbers: 13.60.Le, 14.20.Gk, 25.40-h 

s : 

. The COSY-TOF collaboration [jj has recently measured the w-meson angular distribution in proton-proton colli- 
ng ' sions at an excess energy e = 173 MeV (y/s — 2m p + m u + e) above the uj threshold. Earlier, angular distribution 
of uj- and (fhmesons have been measured by the DISTO collaboration |2j at an energy excess of e = 320 MeV. The 
t-H ■ 0-meson distribution is found to be consistent with a flat distribution, whereas the w-mesons at the both values of e 
k*" I are peaked towards the beam directions. 

The cross sections of vector mesons production enter as an input into transport models for the dilepton emission 
in heavy-ion collisions The direct p(uj) — > e + e~ decay channels give important contributions to the dilepton 

production. The detector efficiency depends thereby on the absolute value and direction of momentum of the dilepton 
pairs. Hence, the angular distribution of the vector mesons produced in nucleon-nucleon collisions affects the counting 
, rate of the dilepton pairs which makes it necessary to keep the angular dependence of the production cross sections 
' theoretically under the control. 

The dilepton spectra in heavy-ion collisions at T = 1 — 5 GeV/ A were measured by the DLS collaboration 0. The 
HADES experiment is presently studying the dilepton spectra in the same energy region in great detail. 

The pion production is the best known inelastic reaction in pp collisions. A model independent partial-wave analysis 
of the experimental data has been performed up to e„ — 20 MeV above the threshold || . The pion angular distribution 
shows an increasing anisotropy with energy. 

The meson production in pp collisions has previously been analyzed using the OBEP models 9] . The most of the 
realistic nucleon-nucleon interactions models have the A(1232) isobar included as a special degree of freedom. It is 
naturallyto assume that with increasing the energy other nucleon resonances also become important. The resonance 
model [H EH 

as a first step treats matrix elements of the higher-mass resonance production as phenomenological 
constants fixed by fitting the available inelastic cross sections. Such an approach is suitable to obtain estimates for 
the importance of various meson production channels and it allows to perform simulations of heavy-ion reactions in 
which the role of nucleon resonances is important. The dominance of nucleon resonances in the intermediate states is 
equivalent to taking nucleon- meson final-state interactions into account [12|. The resonance model fits the p, u), and 
</> total production cross sections rather well 0, 0, 0, 0] . 

The decay rates of nucleon resonances into the dilepton pairs and vector mesons with the same invariant mass 
coincide up to an overall kinematic factor. The A(1232) Dalitz decay is one of the major sources of dileptons in 
heavy- ion collisions at intermediate energies 0, 0, 0- The first correct calculation of that decay was given only 
recently whereas kinematically complete expressions for Dalitz decays of other high-spin resonances were given 
in ^(|. In this paper, we analyze the angular distributions of the uj- and 0-mesons in proton-proton collisions within 
the framework of the resonance model. 

According to the resonance model, the dominant contributions to the meson production cross sections originate from 
reactions pp — > pN* followed by subsequent decays N* — > puj{4>). The first reaction takes place at m p + mjy* < \fs, 
while the second reaction with e.g. an uj production, takes place at m p + m u < tojv». These inequalities can be 
combined to give 1.72 < m N , < 1.89 GeV and 1.72 < m N , < 2.05 GeV, respectively, for the COSY-TOF and DISTO 
experimental conditions. The 0-meson production is sensitive to the mass interval 1.96 < mjy* < 2.04 GeV. 

Either a nucleon resonance falls into these intervals with its pole masses or, if not, it can have a Breit-Wigner leaking 
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TABLE I: j p states of two protons at j > 1 and their 2s+1 lj decompositions satisfying the Pauli principle. 



■p 2s+l, 

2 h 

j + even 1 j j 

j~ odd 3 jj 

j + odd forbidden 

j-even 3 (j - 1), 3 (j + 1) 3 



into the corresponding mass intervals. In the present work we will focus mainly on the first possibility. According to 
Particle Data Group H3 the following resonances can be important: A r *(1720)| + , N*(1900)§ + , Af*(1990)| + , and 

AT* (2000) f . The negative-parity resonances have masses below 1.72 GeV, while A* resonances do not contribute to 
the a; and (f> production due to the isotopic symmetry. 

At e — 173 MeV the resonances A^*(1720)| + and A r *(1900)| + appear at different kinematic conditions: In the 
AT* (1720) decay, the u> is mainly produced at rest. In the pp collisions, the A~*(1720) is produced with a velocity 
v ~ 0.26. The w-meson angular distribution coincides therefore with the angular distribution of the A''* (1720) 
resonance. The maximal angular momentum for the system pN* can be estimated as rj — p* /m w where p* is the cm. 
momentum in the final state, l/(2m,r) is the strong interaction radius and m w is the pion mass. For e = 173 MeV, 
one has r\ ~ 3.4. One can expect that the partial-wave decomposition is truncated at a pN* orbital momentum ~ 3. 

In the A'* (1900) decay the Li is produced with velocity v ~ 0.82, whereas in the pp collisions the AT* (1900) is 
produced mainly at rest. The angular distribution of the aj-mesons coincides therefore with the corresponding uj 
distribution in the Af*(1900) resonance decay. This means that direction of the u> momentum is correlated with 
the resonance spin. The initial protons in the experiments Q, Q are not polarized. The resonance spin is therefore 
correlated with the beam direction. 

Let us first consider the A 7 "* (1900). The angular wave function of the pu) system, appearing due to a J p resonance 
decay, has the form 

*K) ~^A SL Ci£ LIi Cl% ita Y LL .(n ul ). (1) 

SL 2 

The Clebsch-Gordan coefficients C'Jf z „ arc defined with the PDG phase conventions 117b . The partial-wave de- 
composition includes the total spins S — 1/2 and 3/2 of the pto system and the pu orbital momenta L such that 
| J — S\ < L < J + S. Parity conservation gives a selection rule P = (— l) i+1 . The amplitudes Asl for different 
resonances are extracted from the ttA^ inelastic scattering data and/or predicted by the quark models (for a review 
see [l6|'). 

The pA^*(1900) system appears close to the threshold in the S state (77 ~ 0). The initial protons have therefore the 
total angular momentum j = 1 or 2 and positive parity. The state 1 + is forbidden by the Pauli principle, so the only 
possibility is j p = 2 + which is the 1 d,2 state for the initial protons. The higher states are classified in Table 1. 

The transition amplitude pp — ► pN* (1990) can be obtained constructing a rotational scalar out of the initial-state 
spherical harmonics Yjj z (n p ), where n p is the unit vector in direction of the proton beam. The final-state pAf*(1900) 
spin wave function reads: 

M fi~&frjj, Y >j.(*p)- ( 2 ) 

Now, let the quantization axis of the angular momentum be parallel to the proton beam n p . The only non- vanishing 
component of the initial-state wave function is j z = 0. The orbital momentum is perpendicular to the beam and 
its projection to the beam direction is equal to zero. The polarization matrix of a spin- J resonance becomes then 
diagonal. The diagonal elements equal 

PJ.J. = \ C i-j.jj/> ( 3 ) 

and so pj z j z = 1/2 at J z = ±1/2 and, for a higher spin resonance, pj z j z = at \J Z \ > 3/2. Eq.Q is valid for all 
even-j positive-parity initial states. Spin-half resonances are therefore produced unpolarized near threshold, whereas 
higher spin resonances, including the A^*(1900), are polarized. 
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FIG. 1: The angular distribution of the w-mesons in the cm. frame of the colliding protons assuming the reaction goes through 
the N*( 1900)^ resonance. The experimental data from COSY-TOF Q were obtained for an excess energy e = 173 MeV. The 
Pi/2 contribution ~ 1 + 3cos 2 6Z in the pui decay of the AT* (1900) is dominant according to the models MS |T3 |. K [lfl| . and 
eVMD O. The predictions of the models CR [2FJ and SST ^ are shown as well. 



The angular distribution of the cj-mesons produced through the AT* (1900) resonance decays can be found by 
weighting the polarization matrix (J3J with the modulus squared of the puo angular wave function summed over 
the final-state pu spin projections: 



da 
dtt 



E^E 
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The coefficient C 1 " entering the puj wave function (cf. Eq.Q) gives after summation over the n and fi u a decoherent 

sum over the total spin and spin projection of the pu> system. 

The partial- wave amplitudes Als of the TV* (1900) — > pco decay are extracted b y M anley and Saleski from 
ttN multichannel partial- wave analysis and/or quark model predictions by Koniuk |19|. Capstick and Roberts [2fj| . 
and Stassart and Stancu [2l| and/or predicted by the extended Vector Meson Dominance (eVMD) model [T^- The 
resulting angular distributions are shown in FigQ] together with the experimental data 0. 

The models |l6t Il8l Il9l| pr edict a dominant P1/2 wave with an angular dependence ~ 1 + 3cos 2 #* very close to the 
experiment. The models fwL fl8L IT9L l2l | are in very good agreement with the data. Capstick and Roberts [2(| report 
a vanishing p±/2 amplitude and large errors for the p 3 / 2 and f 3 / 2 amplitudes. The dominance of the / 3 / 2 transition 
does not contradict to the model |20| and to the experimental data [l|. This amplitude is shown on Fig^as C-Rf/a/a]- 
The P3/2 transition can apparently be excluded. 

If the iV*(1900) resonance would be unpolarized, it could produce an isotropic w-meson distribution. The distribu- 
tions plotted on FigQ] within the resonance model, are a consequence of the selection rules, the kinematic conditions, 
and the partial- wave content of the iV*(1900) decay amplitude. 

The TV* (1900) width equals T^E = 500 ±80 01 ■ I n general, resonances are produced off-shell with invariant masses 
away from the pole mass. The off-shell A r *(1900)'s move with finite velocities which leads to a smearing of the to 
angular distribution compared to the iV*(1900) decays at rest. To estimate this effect, let us compare velocities of 
the u in the iV*(1900) rest frame and of the iV*(1900) in the c. m. frame of two initial protons. For e = 173 MeV, 
we have = 0.46 and vn* — at m^. = 1.9 GeV and — and «at« = 0.26 at kin* = 1.72 GeV. The ui and N* 
velocities are equal = djj. = 0.23 at tun* — 1-76 GeV, i.e. a 40 MeV above the N* — * pcu decay threshold. The 
approximation in which the TV* (1900) is treated as a particle at rest is therefore reasonable. Hence, we do not expect 
a large smearing effect due to the off-shell production of the iV*(1900). 

Let us now discuss the A^*(1720). The maximum pN* angular momentum at the excess energy e = 173 MeV was 
estimated as ~ 77 = 3.4. Hence the pN* system can appear in orbitally excited states. The partial-wave content of 
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the transition pp — > pN* is unknown. A dominant S-wave component can immediately be excluded since it produces 
an isotropic TV* (1720) and correspondingly also an isotropic u> angular distribution. Let us check if L = 1 and L = 2 
waves are allowed as dominant components. 

There exist eight L = 1 amplitudes: 3 p 3 P , 3 Pi 3 Pi, 3 P2 3 Pu 3 Pi 5 Pi, 3 P2 5 P2, 3 f2 3 P2, 3 h 5 P2, and 3 / 3 5 -P3 with 
quantum numbers 2s+1 lj 2S+1 Lj where s, I, and j are the total spin, orbital momentum, and total angular momentum 
of the initial pp state. S and L are the total spin and orbital momentum of the final pN* state. The angular 
distributions for the isolated transitions can be calculated from 



They are plotted on Fig[2J Since in the resonance rest frame the w-meson is produced at rest, the u> momentum has 
in the pp cm. frame the same direction as the iV*(1720) momentum. Thus we have replaced in (JJJ njv* with n w . 

The wave Nos. 6, 8, 3, and 2 with quantum numbers 3 fz 3 P^, 3 f3 5 P3, 3 P2 3 Pi, and 3 pi 3 Pi resemble qualitatively 
the measured distribution, but yet no conclusions can be drawn on the actual importance of these waves. The ratio 
between maximal and minimal values of the differential cross section for those waves is below a factor of 3, whereas 
the experimental ratio is close to a factor of 4. The resonance A* (1900) alone gives a stronger angular dependence, 
being in slightly better agreement with the data. 

There exist four L = 2 transition amplitudes: 1 d-2 3 Di, 1 sq 5 Do, 1 d,2 5 D2, and 1 g^D±. The corresponding w-meson 
angular distributions are shown in Fig|5j The 1 g<$ 5 D± distribution resembles the experimental one, but has a plateau 
at cos6^ ~ 0. The pattern of other waves is quite different from the observed one. 

The transitions L = and L = 1 Nos. 1, 4, 5, and 7 and L = 2 Nos. 1, 2, and 3 can apparently be excluded as 
dominant ones. The other four P-waves and the one D-wave transition can be large. These conclusions, however, are 
not stable against effects from the finite A* (1720) width of T^l = 150 ± 50 MeV 

The off-shell iV*(1720) production results into an to distribution representing a convolution of the 7V*(1720) and 
oj angular distributions. As we estimated, at m^* = 1-76 GeV the resonance and u> velocities are equal, so the 
A*(1720) tail m^* > 1-76 GeV generates roughly the A*(1720) — ■> Nuo decay distribution, whereas the invariant 
masses 1.72 < tun* < 1-76 GeV correspond to one of the distributions plotted on Fig|3] The domain of 40 MeV 
above the u> threshold is not lar ge c ompared to the resonance width. Hence the pattern of the decay A* (1720) — * Nu 
is important. According to [Ta Il9| the p\/2 wave ~ 1 + 3cos 2 #* of the A* (1720) — > Nuj decay is dominant. Other 
authors give negligible uj couplings for the A* (1720). 

The A*(1720) production cross section is parameterized as an S'-wave process This resonance provides 10 — 20% 
of the total ppu> cross section [l3| . Such a value does not contradict to the cu- meson angular distribution even in the 
sharp resonance limit. The angular distribution data do not provide stringent constraints to the S'-wave part of the 
V*(1720) production. 

The A* (1900) is reliably predicted by quark models, but has currently a lower experimental status (**) than the 
A*(1720). The A*(1900) cross section is unknown and its contribution to the vector mesons production is set equal 
to zero 0,E3>E1- In the region y/s ~ 2.8 GeV the vector meson production, according to 0^E10> is determined 
by tails of resonances of smaller masses and by the A* (1720). If we attribute the total cross section to the TV* (1900) 
alone, we get a(pp -> pA*(1900))5(A*(1900) -> pu>) ~ 30 fib. Using the estimate of [3 T N * pLJ ~ 60 MeV and the 
PDG value of Tjfit ~ 500 MeV pj, one gets a(pp -> pA*(1900)) - 0.3 mb at y/s ~ 2.8 GeV. This is below the s-wave 
unitarity limit of 4ir/k 2 ~ 4 mb. 

Under the conditions of the DISTO experiment, we have = and ujv* = 0.35 at mjv» = 1.72 GeV and 
v u = 0.56 and vn* = at mjy* = 2.05 GeV. The resonance and u> velocities at e = 320 MeV are equal to each other 
Vu> = vn* = 0.3 at mjv» = 1.8 GeV. 

The approximation v u ~ is good for the resonance A*(1720). The angular distribution due to the A*(1720) 
production is the same as in Fig|2] The L — 2 wave 1 54 5 Z?4 is the only wave that fits the DISTO w-meson data. 

The approximation vat. ~ is good for the resonances A* (1900), A* (1990), and A* (2000). The angular distribution 
due to the A* (1900) production is the same as in Fig^ The pattern of the distribution is correct, but the A* (1900) 
does not develop a plateau at cos#* ~ which is required by the DISTO data. The angular distributions of the 
V*(1990) and A*(2000) are shown on Fig|31 The model [lij gives results close to the observations in both cases, the 
other models give reasonable angular dependences in the case of a A* (1990) dominance but stand in contradiction to 
a possible A* (2000) dominance. 

The 0-meson angular distribution at = 2.98 GeV, i.e. 85 MeV above the 0-meson production threshold, is 
found to be consistent with an isotropic distribution 0- The maximal cf> momentum can be estimated in the cm. 
frame to be p™ ax = 0.33 GeV and the maximal <fi orbital momentum relative to the pp pair ~ r/^ = p™ ax /m w = 2.4. 
The S-wave (/>-meson state allows j p = 1~ pp<f> final states, whereas the higher orbital states are apparently excluded 
by the data which is somewhat surprising. 
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T resonance in transitions 2a+1 lj 2S+1 Lj where 



FIG. 2: The normalized L = 1 and L = 2 angular distributions of the TV* (1720)^ 
s, I, and j are the total spin, orbital momentum, and total angular momentum in the initial pp state. S and L are the total 
spin and orbital momentum in the final pN* state. The numbers 1 to 8, attached to the L = 1 curves, correspond to the waves 
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o, 3 pi 3 Pi, 3 p 2 3 Pi, 3 Pi 5 Pi, 3 p 2 5 P 2 , 3 f 2 3 P 2 , 3 /2 5 ^2, and 3 / 3 5 P 3 , respectively. The numbers 1 to 4, attached to the L = 2 
curves, correspond to the waves 1 d 2 3 D 2 , 1 s 5 D , 1 d 2 5 D 2 , and 1 g^D i . 
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FIG. 3: The angular distribution of the oj-mesons in the cm. frame of the colliding protons assuming that the reaction proceeds 
through TV* (1990) | + and TV* (2000) | + decays. The models K 19], CR 20"], SST 2Xj, and eVMD are compared to the 
DISTO data at Vs = 2.98 GeV @] 



The u — <f> mixing is sufficient to reproduces the c6-meson production cross section at </s = 2.98 GeV Such a 
mechanism, owing to the mass dependence of the NN*lo coupling constants, yields in the nucleon resonance decays 
identical angular distributions for <f>- and w-mesons. To construct a T final state out of the proton and a resonance 
in a S'-wave, one needs negative-parity resonances. The well established negative-parity resonances iV*(1535)i , 

TV* (1650) i~, and 7V*(1700)f" with non-vanishing couplings to w-mesons have low masses and can only contribute 
to the reaction through their Breit-Wigner tails. These resonances have dominant si/ 2 branchings [la. fl9l| and thus 

their tails would produce an isotropic cross section. The resonances 7V*(2090)| and A r *(2080)| are listed by PDG 
[l7| with ratings (*) and (**), respectively. Quark models predict several negative-parity states with masses around 
2000 MeV but according to ref. 20] the si/ 2 decays of these resonances are subdominant. 

The J = \ resonances in the j p = 1~ state give isotropic distributions. In the reference frame where the quantization 
axis of the angular momentum is parallel to the proton beam, the polarization matrix of the spin- J resonance can be 
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FIG. 4: The angular distribution of the 0-mesons in the cm. frame of the colliding protons measured by the DISTO collabora- 
tion at t/s = 2.98 GeV Q is compared to the flat distribution (solid line) generated by the dominant N* (2090) 5 contribution 
in the reaction mechanism. 



found to be 

p^ = ^^\ c ^: M J- (6) 

For s — 0, we recover Eq.(3). Eq.@ is valid for all 2s+1 lj states of two protons. The spin-half nucleon resonances 
are produced unpolarized. The diagonal elements of the J — 3/2 polarization matrix equal pj r j z =1/8 and 3/8 for 
J z = ±1/2 and ±3/2, respectively. So, the -/V*(2090)^ dominance in the pp — > pp<f> or the dominance of any other 
spin-half nucleon resonance would result in a flat 0-meson distribution. The high-spin resonances are polarized and 
their dominance can apparently be excluded (except for J = - resonances with the dominant Sx/2 decay modes). 

Since the energy released in the process pp — > pp<f> is small, it is instructive to compare the <f> production with the 
near-threshold pion production which is well studied at < 20 MeV @. The value of 77^ = p™ ax /m, is small and 
varies in the limits 0.07 — 0.5. In the limit T] v — > 0, only the S-wave survives, so the distribution must be symmetric. 
The large S-wave of the NNir system e.g. at r) v ~ 0.22 implies an almost symmetric angular distribution which 
is reported 8]. Within the resonance model, the dominance of the S-wave implies apparently the dominance near 
the threshold of tails of the negative-parity nucleon resonances produced in an S'-wave together with the nucleon. 
The importance of heavy resonances near the pion threshold is emphasized in 22] . With increasing t] n , the pion 
distribution develops an anisotropy that becomes essential at rj n ~ 0.5. It can be attributed at least partially to 
the A(1232) which appears in the S'-wave with another nucleon and generates the pion distribution ~ 1 + 3cos 2 #*, 
the same as the w-meson distribution in the P1/2 wave with the proton, since the A(1232) has the same P1/2 decay 
mode and the same spin as the jV*(1900) resonance. Note that Eq.(@J has no dependence on the constituent spins, 
only the total spin S is important. At = 20 MeV, the ratio between the j p = 1~ wave cross section that gives 
the flat distribution and the j p — 2 + wave cross section that gives the ~ 1 + 3cos 2 9* distribution can be estimated 
from the angular distribution alone to give ~ 3 : 2 which is almost two times lower than the ratio of ~ 5 : 2 given 
by 8J. There is therefore space for contributions due to the asymmetry of other resonances and/or for higher orbital 
states of the iVA(1232) system 23] . In the pp — > ppcf> reaction, while the 0-meson is nonrelativistic, the maximal 
value of the orbital momentum ~ rj^ — p , ( p ax /m n is several times larger than in the pion production studied in 0. 
These two near-threshold processes appear to be different with respect to the convergence rate of the partial-wave 
expansion. The isotropic distribution observed by the DISTO collaboration can be interpreted as an evidence for a 

strong dominance of the TV* (2090) ^ or any other spin-half negative-parity resonance. We would otherwise expect on 
the pure kinematic grounds an isotropic angular distribution starting from the lo production threshold of ^/s^ = 2.90 
GeV up to 2.92 GeV where rj^ < 1 and an increased anisotropy beyond 2.92 GeV where higher orbital waves should 
come into play. 

The angular differential cross sections provide sensitive tests for the reaction mechanisms. The present data do 
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not offer enough information for a complete partial-wave analysis. The resonance model which we discussed can be 
considered as first step to construct a microscopic NN — ► NN* transition potential usin g th e i-channel exchange 
OBEP models or proposed recently an s-channel dibaryon exchange NN interaction model [2J] . 

The Y*(1720)f and V*(1900)§ resonances both have the dominant p 1 / 2 decay Nlo modes. If velocities of these 
resonances in the c. m. frame of two protons are not high compared to the lo velocities in the N* rest frames, which is 
a good approximation also for a significant off-shell spectral part of the TV* (1720), the N* — ► Nlo decays produce an 
angular distribution very close to the COSY-TOF data. The mechanism involving these two resonances could be the 
dominant one in the lo production at ^/s = 2.83 GeV. The several P and the one D wave of the W*(1720) system 
generate also lo patterns close to the experimentally observed one. 

At ^fs = 2.98 GeV a dominance of TV* (1990) f + and also TV* (2000) | + within the model HJ gives a reasonable 
description of the DISTO angular distribution data. The resonance V*(1900)| + at rest generates still the right 
pattern but without a plateau at cos^* ~ 0. The TV* (1720) generates the right pattern through the 1 gi 5 D^ transition. 

The (/j-meson data can be explained by the dominance of the iV*(2090)i . It is not clear, however, why higher 
orbital states of the <fi and, respectively, high-spin nucleon resonances which give anisotropic distributions should not 
be involved in the </>-meson production. 

We formulated also a general theorem according to which from every 2s+1 lj initial state of two unpolarized protons, 
spin-half nucleon resonances are produced at the threshold being unpolarized, whereas high-spin nucleon resonances 
are produced being polarized. The absence of any structure in the meson angular distribution indicates the possible 
dominance of a spin-half nucleon resonance, whereas the presence of a structure indicates the possible dominance of 
a J > 3/2 nucleon resonance. 

A complete partial wave analysis requires angular distributions from experiments with polarized beams. At the 
moment several competitive mechanisms can explain the w-meson production data while the (/)-meson distribution 
cannot be interpreted so easily. 

The authors are grateful to V. I. Kukulin for reading the manuscript and valuable remarks. M.I.K. and B.V.M. 
wish to acknowledge kind hospitality at the University of Tuebingen. This work is supported by DFG grant No. 436 
RUS 113/721/0-1 and RFBR grant No.03-02-04004. 
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